-----------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\jgw12\Dropbox\Research\PersPartyBook\Data\FKTW-reproduction\C
> h5.log
  log type:  text
 opened on:  26 Jul 2023, 17:03:03

. 
.         ******************************
.         **** Set directory, seed *****
.         ******************************
.                 set more off 

.                 set matsize 1000
set matsize ignored.
    Matrix sizes are no longer limited by c(matsize) in modern Statas.  Matrix
    sizes are now limited by edition of Stata.  See limits for more details.

.                 global seed ="984353"

.                 set scheme plotplain

.                 cd "$dir"
C:\Users\jgw12\Dropbox\Research\PersPartyBook\Data\FKTW-reproduction

. 
.         ********************************
.         **** Load data, vars trans *****
.         ********************************
.                 use pers-use,clear

.                 qui reg v2x_polyarchy  create ld if persparty~=.

.                 gen sample = e(sample)==1

.                 egen c=count(year) if e(sample)==1,by(lid) 

.          
.                 * Variables and transformations *
.                 gen election =  v2xel_elecparl==1 | v2xel_elecpres==1

.                 gen seat50 = v2paseats>49 if v2paseats~=.
(174 missing values generated)

.                 gen time = year-1990

.                 
.                 * Standardize outcome variables *
.                 local var = "vburcap v2lgotovst v2lginvstp v2juhcind v2xlg_legcon
>  v2x_jucon"

.                 foreach v of local var {
  2.                         qui sum `v'
  3.                         qui replace `v'=`v'-r(mean)
  4.                         qui sum `v'
  5.                         qui replace `v'=`v'/r(sd)
  6.                 }

.  
.                 sum vburcap v2lgotovst v2lginvstp v2juhcind v2xlg_legcon v2x_juco
> n

    Variable |        Obs        Mean    Std. dev.       Min        Max
-------------+---------------------------------------------------------
     vburcap |      2,392   -3.37e-09           1  -2.072982   1.991608
  v2lgotovst |      2,384    2.65e-17           1  -3.237818   1.958456
  v2lginvstp |      2,385    1.07e-16           1   -3.41118   2.153816
   v2juhcind |      2,392   -7.74e-17           1  -3.047219   2.102983
v2xlg_legcon |      2,384    6.52e-17           1  -3.540382   1.090339
-------------+---------------------------------------------------------
   v2x_jucon |      2,392   -1.14e-17           1  -3.467651   1.101739

.                 * Initial levels * Note year before election that brings leader t
> o power
.                 local var = "polar v2xps_party vburcap v2lgotovst v2lginvstp v2xl
> g_legcon v2x_jucon v2juhcind "

.                 foreach v of local var {
  2.                         qui tsset cowcode year
  3.                         qui gen o`v'=(l2`v') if minyr==year
  4.                         qui egen i`v'=max(o`v'),by(lid)
  5.                         qui gen d`v'=d.`v'
  6.                 }

.                 
.                 * Globals *
.                 global d="persparty"

.             global x1="ld ivdem time"           

.                 
.                 * Strong autocorrelation decreases with differenced DV *
.                 local var = "v2xlg_legcon v2lginvstp v2x_jucon v2juhcind vburcap 
> v2lgotovst"

.                 foreach v of local var {
  2.                            xtset cowcode year
  3.                            qui prais `v' persparty $x1 i`v' $t, cluster(lid)
  4.                                         di `e(rho)'
  5.                            qui prais `v' persparty $x1 l1`v' $t, cluster(lid)
  6.                                         di `e(rho)'
  7.                            qui prais `v' persparty $x1 l1`v' l2`v' $t, cluster
> (lid)
  8.                                         di `e(rho)'
  9.                  }

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit
.98956659
.12627743
.01077828

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit
.9847208
.01396222
-.04634113

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit
1.0004904
.08405653
-.00893034

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit
.98234223
.02781819
-.01246123

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit
.99511672
.06604835
.00226108

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit
.98647532
-.03589383
-.07452498

.                  
.                  **************************************
.                  ********* Reported estimates *********
.                  **************************************
.         * Legislative constraint by party creation (Table 5.1) *
.                  tab seat50

     seat50 |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,601       72.18       72.18
          1 |        617       27.82      100.00
------------+-----------------------------------
      Total |      2,218      100.00

.                  reg v2xlg_legcon create,cluster(lid)

Linear regression                               Number of obs     =      2,384
                                                F(1, 591)         =      13.56
                                                Prob > F          =     0.0003
                                                R-squared         =     0.0347
                                                Root MSE          =     .98272

                                  (Std. err. adjusted for 592 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
      create |  -.4115543   .1117695    -3.68   0.000    -.6310681   -.1920404
       _cons |   .1180802   .0585715     2.02   0.044     .0030465    .2331138
------------------------------------------------------------------------------

.                  est store l1

.                  reghdfe v2xlg_legcon create,a(cowcode)cluster(lid)
(dropped 4 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,380
Absorbing 1 HDFE group                            F(   1,    587) =       5.74
Statistics robust to heteroskedasticity           Prob > F        =     0.0169
                                                  R-squared       =     0.8498
                                                  Adj R-squared   =     0.8430
                                                  Within R-sq.    =     0.0148
Number of clusters (lid)     =        588         Root MSE        =     0.3947

                                  (Std. err. adjusted for 588 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
      create |   -.140613   .0586899    -2.40   0.017    -.2558808   -.0253453
       _cons |   .0441127   .0189624     2.33   0.020     .0068704     .081355
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       102           0         102     |
-----------------------------------------------------+

.                  est store l2

.                  reghdfe v2xlg_legcon create seat50,a(cowcode)cluster(lid)
(dropped 4 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,206
Absorbing 1 HDFE group                            F(   2,    537) =      11.85
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.8633
                                                  Adj R-squared   =     0.8568
                                                  Within R-sq.    =     0.0739
Number of clusters (lid)     =        538         Root MSE        =     0.3747

                                  (Std. err. adjusted for 538 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
      create |   -.146423   .0585838    -2.50   0.013    -.2615045   -.0313415
      seat50 |  -.2660372   .0549711    -4.84   0.000    -.3740219   -.1580524
       _cons |   .1488084   .0253695     5.87   0.000     .0989728     .198644
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       100           0         100     |
-----------------------------------------------------+

.                  est store l3

.                  xi:reghdfe v2xlg_legcon  create if seat50==1,a(cowcode)cluster(l
> id)
(dropped 7 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =        609
Absorbing 1 HDFE group                            F(   1,    151) =      10.27
Statistics robust to heteroskedasticity           Prob > F        =     0.0017
                                                  R-squared       =     0.8847
                                                  Adj R-squared   =     0.8728
                                                  Within R-sq.    =     0.0747
Number of clusters (lid)     =        152         Root MSE        =     0.4116

                                  (Std. err. adjusted for 152 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
      create |  -.5365399   .1674305    -3.20   0.002    -.8673489   -.2057309
       _cons |  -.2458075   .0517643    -4.75   0.000    -.3480832   -.1435317
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        57           0          57     |
-----------------------------------------------------+

.                  est store l4

.                  xi:reghdfe v2xlg_legcon  create if seat50==0,a(cowcode)cluster(l
> id)
(dropped 6 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      1,588
Absorbing 1 HDFE group                            F(   1,    452) =       0.80
Statistics robust to heteroskedasticity           Prob > F        =     0.3709
                                                  R-squared       =     0.9056
                                                  Adj R-squared   =     0.8996
                                                  Within R-sq.    =     0.0021
Number of clusters (lid)     =        453         Root MSE        =     0.2758

                                  (Std. err. adjusted for 453 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
      create |  -.0369776   .0412821    -0.90   0.371    -.1181062     .044151
       _cons |   .2041807   .0142492    14.33   0.000     .1761778    .2321835
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        94           0          94     |
-----------------------------------------------------+

.                  est store l5

.                 
.         * Legislative constraint on executive *
.                  reg v2xlg_legcon l1v2xlg_legcon l2v2xlg_legcon $d $x1 seat50,clu
> ster(lid)

Linear regression                               Number of obs     =      2,193
                                                F(7, 537)         =    2669.88
                                                Prob > F          =     0.0000
                                                R-squared         =     0.9499
                                                Root MSE          =     .22223

                                    (Std. err. adjusted for 538 clusters in lid)
--------------------------------------------------------------------------------
               |               Robust
  v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
---------------+----------------------------------------------------------------
l1v2xlg_legcon |   5.270729   .1941339    27.15   0.000     4.889374    5.652084
l2v2xlg_legcon |  -.5555872   .1845327    -3.01   0.003    -.9180817   -.1930926
     persparty |  -.0884558   .0290251    -3.05   0.002    -.1454725   -.0314391
            ld |  -.0186393    .007379    -2.53   0.012    -.0331346    -.004144
         ivdem |   .1370543   .0532681     2.57   0.010     .0324149    .2416938
          time |  -.0007767   .0005822    -1.33   0.183    -.0019203    .0003669
        seat50 |  -.0394189   .0129583    -3.04   0.002    -.0648741   -.0139637
         _cons |  -3.579826   .0471702   -75.89   0.000    -3.672487   -3.487166
--------------------------------------------------------------------------------

.                  est store c1

.                  reghdfe v2xlg_legcon l1v2xlg_legcon l2v2xlg_legcon $d $x1 seat50
> ,a(cowcode)cluster(lid)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,190
Absorbing 1 HDFE group                            F(   7,    534) =     133.37
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9571
                                                  Adj R-squared   =     0.9549
                                                  Within R-sq.    =     0.7067
Number of clusters (lid)     =        535         Root MSE        =     0.2099

                                    (Std. err. adjusted for 535 clusters in lid)
--------------------------------------------------------------------------------
               |               Robust
  v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
---------------+----------------------------------------------------------------
l1v2xlg_legcon |   4.449437   .3114737    14.29   0.000     3.837573    5.061301
l2v2xlg_legcon |  -.5134696   .1837635    -2.79   0.005    -.8744576   -.1524815
     persparty |  -.1631885   .0425311    -3.84   0.000    -.2467373   -.0796396
            ld |  -.0517201   .0200164    -2.58   0.010    -.0910406   -.0123996
         ivdem |   .4486176   .1732848     2.59   0.010      .108214    .7890211
          time |  -.0002368   .0009299    -0.25   0.799    -.0020634    .0015898
        seat50 |  -.0621665   .0194035    -3.20   0.001    -.1002831     -.02405
         _cons |  -3.064583   .1440347   -21.28   0.000    -3.347527   -2.781639
--------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       100           0         100     |
-----------------------------------------------------+

.                  est store d1

.                  xi:reghdfe v2xlg_legcon l1v2xlg_legcon l2v2xlg_legcon i.seat50*$
> d $x1,a(cowcode)cluster(lid)
i.seat50          _Iseat50_0-1        (naturally coded; _Iseat50_0 omitted)
i.seat50*pers~y   _IseaXpersp_#       (coded as above)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,190
Absorbing 1 HDFE group                            F(   8,    534) =     119.74
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9571
                                                  Adj R-squared   =     0.9549
                                                  Within R-sq.    =     0.7070
Number of clusters (lid)     =        535         Root MSE        =     0.2099

                                    (Std. err. adjusted for 535 clusters in lid)
--------------------------------------------------------------------------------
               |               Robust
  v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
---------------+----------------------------------------------------------------
l1v2xlg_legcon |   4.435022   .3141316    14.12   0.000     3.817937    5.052107
l2v2xlg_legcon |  -.5123236   .1836207    -2.79   0.005    -.8730311   -.1516162
    _Iseat50_1 |  -.0175852     .03488    -0.50   0.614    -.0861041    .0509336
     persparty |  -.1360514   .0417791    -3.26   0.001     -.218123   -.0539799
 _IseaXpersp_1 |  -.0866842   .0859526    -1.01   0.314    -.2555309    .0821625
            ld |  -.0534973   .0202753    -2.64   0.009    -.0933263   -.0136682
         ivdem |   .4535983   .1722808     2.63   0.009     .1151671    .7920295
          time |  -.0000888    .000963    -0.09   0.927    -.0019805     .001803
         _cons |  -3.068773   .1427989   -21.49   0.000     -3.34929   -2.788257
--------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       100           0         100     |
-----------------------------------------------------+

.                  est store s1

.                  interflex v2xlg_legcon $d v2paseat l1v2xlg_legcon l2v2xlg_legcon
>  $x1,fe(cowcode)cluster(lid)nbin(4)
p value of Wald test: 0.0009

.                  mat list r(estBin)

r(estBin)[4,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1       17.15  -.10144549   .05507225  -.20938513   .00649414
r2        32.3  -.18292994    .0559344  -.29255936  -.07330052
r3        45.1  -.15807934   .05574315  -.26733391  -.04882477
r4        56.6  -.24967134   .08934866   -.4247915  -.07455118

.                  reghdfe v2xlg_legcon l1v2xlg_legcon l2v2xlg_legcon $d $x1 seat50
> ,a(cowcode year)cluster(lid)
(dropped 3 singleton observations)
note: time is probably collinear with the fixed effects (all partialled-out values 
> are close to zero; tol = 1.0e-09)
(MWFE estimator converged in 6 iterations)
note: time omitted because of collinearity

HDFE Linear regression                            Number of obs   =      2,190
Absorbing 2 HDFE groups                           F(   6,    534) =     155.95
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9577
                                                  Adj R-squared   =     0.9549
                                                  Within R-sq.    =     0.7090
Number of clusters (lid)     =        535         Root MSE        =     0.2098

                                    (Std. err. adjusted for 535 clusters in lid)
--------------------------------------------------------------------------------
               |               Robust
  v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
---------------+----------------------------------------------------------------
l1v2xlg_legcon |   4.494107   .3058309    14.69   0.000     3.893328    5.094886
l2v2xlg_legcon |  -.5499813   .1823859    -3.02   0.003    -.9082631   -.1916995
     persparty |  -.1644322   .0428957    -3.83   0.000    -.2486972   -.0801672
            ld |  -.0503229   .0200718    -2.51   0.012    -.0897521   -.0108936
         ivdem |    .445844   .1670108     2.67   0.008     .1177653    .7739227
          time |          0  (omitted)
        seat50 |  -.0601879   .0190915    -3.15   0.002    -.0976916   -.0226843
         _cons |  -3.077236   .1467283   -20.97   0.000    -3.365472   -2.789001
--------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       100           0         100     |
        year |        30           1          29     |
-----------------------------------------------------+

.                  est store f1

.                  * Not different in pres v parl *
.                  xi:reghdfe v2xlg_legcon l1v2xlg_legcon l2v2xlg_legcon i.pres*$d 
> $x1,a(cowcode)cluster(lid)
i.pres            _Ipres_0-1          (naturally coded; _Ipres_0 omitted)
i.pres*perspa~y   _IpreXpersp_#       (coded as above)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,364
Absorbing 1 HDFE group                            F(   8,    584) =     116.30
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9558
                                                  Adj R-squared   =     0.9537
                                                  Within R-sq.    =     0.7076
Number of clusters (lid)     =        585         Root MSE        =     0.2141

                                    (Std. err. adjusted for 585 clusters in lid)
--------------------------------------------------------------------------------
               |               Robust
  v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
---------------+----------------------------------------------------------------
l1v2xlg_legcon |   4.537996   .2752066    16.49   0.000     3.997481    5.078511
l2v2xlg_legcon |   -.535868   .1675804    -3.20   0.001    -.8650017   -.2067343
      _Ipres_1 |  -.0387061   .0544296    -0.71   0.477    -.1456076    .0681955
     persparty |   -.153869   .0585827    -2.63   0.009    -.2689274   -.0388106
 _IpreXpersp_1 |   .0082277     .07573     0.11   0.914    -.1405087    .1569641
            ld |  -.0593323   .0193752    -3.06   0.002    -.0973859   -.0212787
         ivdem |   .4846174   .1659776     2.92   0.004     .1586317    .8106031
          time |    .000083    .000968     0.09   0.932    -.0018182    .0019842
         _cons |  -3.134178   .1368346   -22.90   0.000    -3.402926    -2.86543
--------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       102           0         102     |
-----------------------------------------------------+

.                  xi:reghdfe v2xlg_legcon l1v2xlg_legcon l2v2xlg_legcon $d $x1 if 
> pres==1,a(cowcode)cluster(lid)
(dropped 2 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      1,182
Absorbing 1 HDFE group                            F(   6,    248) =     102.54
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9549
                                                  Adj R-squared   =     0.9524
                                                  Within R-sq.    =     0.7524
Number of clusters (lid)     =        249         Root MSE        =     0.2339

                                    (Std. err. adjusted for 249 clusters in lid)
--------------------------------------------------------------------------------
               |               Robust
  v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
---------------+----------------------------------------------------------------
l1v2xlg_legcon |    4.65604   .3558194    13.09   0.000     3.955226    5.356853
l2v2xlg_legcon |  -.5072163   .2162514    -2.35   0.020    -.9331399   -.0812928
     persparty |  -.1362729   .0529209    -2.58   0.011    -.2405047   -.0320411
            ld |  -.0745599   .0275013    -2.71   0.007    -.1287257   -.0203941
         ivdem |   .3676018   .2570352     1.43   0.154    -.1386484     .873852
          time |   .0006796   .0016516     0.41   0.681    -.0025733    .0039325
         _cons |  -3.204422   .2001653   -16.01   0.000    -3.598663   -2.810181
--------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        57           0          57     |
-----------------------------------------------------+

.                  xi:reghdfe v2xlg_legcon l1v2xlg_legcon l2v2xlg_legcon $d $x1 if 
> pres==0,a(cowcode)cluster(lid)
(dropped 1 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      1,182
Absorbing 1 HDFE group                            F(   6,    337) =      55.55
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9437
                                                  Adj R-squared   =     0.9410
                                                  Within R-sq.    =     0.5990
Number of clusters (lid)     =        338         Root MSE        =     0.1904

                                    (Std. err. adjusted for 338 clusters in lid)
--------------------------------------------------------------------------------
               |               Robust
  v2xlg_legcon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
---------------+----------------------------------------------------------------
l1v2xlg_legcon |   3.961599   .3992261     9.92   0.000      3.17631    4.746888
l2v2xlg_legcon |  -.5640339   .2491253    -2.26   0.024     -1.05407   -.0739974
     persparty |  -.1756725   .0604795    -2.90   0.004    -.2946374   -.0567075
            ld |  -.0348953   .0263021    -1.33   0.186    -.0866324    .0168418
         ivdem |    .786173   .2793871     2.81   0.005     .2366106    1.335735
          time |   .0000491   .0010145     0.05   0.961    -.0019464    .0020447
         _cons |  -2.899426   .2094767   -13.84   0.000    -3.311473    -2.48738
--------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        49           0          49     |
-----------------------------------------------------+

.                  
.         * Judicial constraint *
.                  reg v2x_jucon l1v2x_jucon l2v2x_jucon $d $x1,cluster(lid)

Linear regression                               Number of obs     =      2,386
                                                F(6, 589)         =    7361.06
                                                Prob > F          =     0.0000
                                                R-squared         =     0.9689
                                                Root MSE          =     .17659

                                  (Std. err. adjusted for 590 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
   v2x_jucon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
 l1v2x_jucon |   5.060097   .1528339    33.11   0.000     4.759931    5.360262
 l2v2x_jucon |    -.39025    .159439    -2.45   0.015    -.7033882   -.0771118
   persparty |  -.0657903   .0250465    -2.63   0.009    -.1149816    -.016599
          ld |  -.0104834   .0056246    -1.86   0.063    -.0215302    .0005633
       ivdem |   .1331898   .0475439     2.80   0.005     .0398135     .226566
        time |  -.0003153   .0004873    -0.65   0.518    -.0012722    .0006417
       _cons |   -3.58676   .0324968  -110.37   0.000    -3.650584   -3.522936
------------------------------------------------------------------------------

.                  est store c2

.                  reghdfe v2x_jucon l1v2x_jucon l2v2x_jucon $d $x1,a(cowcode)clust
> er(lid)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,383
Absorbing 1 HDFE group                            F(   6,    586) =     211.91
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9730
                                                  Adj R-squared   =     0.9717
                                                  Within R-sq.    =     0.7215
Number of clusters (lid)     =        587         Root MSE        =     0.1678

                                  (Std. err. adjusted for 587 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
   v2x_jucon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
 l1v2x_jucon |   4.502357   .1630625    27.61   0.000     4.182099    4.822615
 l2v2x_jucon |  -.4040856    .146643    -2.76   0.006    -.6920954   -.1160757
   persparty |  -.0775784   .0310338    -2.50   0.013    -.1385294   -.0166275
          ld |  -.0438107   .0158808    -2.76   0.006    -.0750009   -.0126206
       ivdem |   .2743892   .1619247     1.69   0.091    -.0436343    .5924127
        time |   .0003841   .0008451     0.45   0.650    -.0012757    .0020439
       _cons |  -3.155178   .1152239   -27.38   0.000     -3.38148   -2.928876
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       103           0         103     |
-----------------------------------------------------+

.                  est store d2

.                  xi:reghdfe v2x_jucon l1v2x_jucon l2v2x_jucon i.seat50*$d $x1,a(c
> owcode)cluster(lid)
i.seat50          _Iseat50_0-1        (naturally coded; _Iseat50_0 omitted)
i.seat50*pers~y   _IseaXpersp_#       (coded as above)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,209
Absorbing 1 HDFE group                            F(   8,    536) =     166.77
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9729
                                                  Adj R-squared   =     0.9716
                                                  Within R-sq.    =     0.7313
Number of clusters (lid)     =        537         Root MSE        =     0.1643

                                   (Std. err. adjusted for 537 clusters in lid)
-------------------------------------------------------------------------------
              |               Robust
    v2x_jucon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
--------------+----------------------------------------------------------------
  l1v2x_jucon |   4.426162   .1640919    26.97   0.000      4.10382    4.748504
  l2v2x_jucon |  -.3423836   .1601522    -2.14   0.033    -.6569865   -.0277808
   _Iseat50_1 |   .0388159   .0310816     1.25   0.212    -.0222408    .0998727
    persparty |  -.0279127   .0312916    -0.89   0.373    -.0893819    .0335565
_IseaXpersp_1 |  -.1749908   .0629068    -2.78   0.006    -.2985649   -.0514167
           ld |  -.0419608    .017621    -2.38   0.018    -.0765755    -.007346
        ivdem |   .2148214   .1583962     1.36   0.176    -.0963321    .5259749
         time |   .0005779   .0009225     0.63   0.531    -.0012343    .0023901
        _cons |  -3.114475    .115433   -26.98   0.000    -3.341231   -2.887718
-------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       101           0         101     |
-----------------------------------------------------+

.                  est store s5

.                  interflex v2x_jucon $d v2paseat l1v2x_jucon l2v2x_jucon $x1,fe(c
> owcode year)cluster(lid)nbin(4)
p value of Wald test: 0.0318

.                  mat list r(estBin)

r(estBin)[4,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        17.2  -.03853588   .03842949  -.11385629   .03678454
r2        32.5  -.00413273   .03998916  -.08251004   .07424458
r3        45.1  -.06197932   .03965132  -.13969448   .01573585
r4       56.65  -.25481343   .06255801  -.37742488  -.13220198

.                  reghdfe v2x_jucon l1v2x_jucon l2v2x_jucon $d $x1,a(cowcode)clust
> er(lid)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,383
Absorbing 1 HDFE group                            F(   6,    586) =     211.91
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9730
                                                  Adj R-squared   =     0.9717
                                                  Within R-sq.    =     0.7215
Number of clusters (lid)     =        587         Root MSE        =     0.1678

                                  (Std. err. adjusted for 587 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
   v2x_jucon | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
 l1v2x_jucon |   4.502357   .1630625    27.61   0.000     4.182099    4.822615
 l2v2x_jucon |  -.4040856    .146643    -2.76   0.006    -.6920954   -.1160757
   persparty |  -.0775784   .0310338    -2.50   0.013    -.1385294   -.0166275
          ld |  -.0438107   .0158808    -2.76   0.006    -.0750009   -.0126206
       ivdem |   .2743892   .1619247     1.69   0.091    -.0436343    .5924127
        time |   .0003841   .0008451     0.45   0.650    -.0012757    .0020439
       _cons |  -3.155178   .1152239   -27.38   0.000     -3.38148   -2.928876
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       103           0         103     |
-----------------------------------------------------+

.                  est store f2

.  
.         * Executive oversight * 
.                  reg v2lgotovst l1v2lgotovst l2v2lgotovst $d $x1,cluster(lid)

Linear regression                               Number of obs     =      2,367
                                                F(6, 587)         =    4260.61
                                                Prob > F          =     0.0000
                                                R-squared         =     0.9480
                                                Root MSE          =     .22786

                                  (Std. err. adjusted for 588 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
  v2lgotovst | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
l1v2lgotovst |   .9246613   .0453492    20.39   0.000     .8355948    1.013728
l2v2lgotovst |   .0003187   .0443272     0.01   0.994    -.0867406     .087378
   persparty |  -.0845375   .0285075    -2.97   0.003    -.1405266   -.0285483
          ld |  -.0220748     .00733    -3.01   0.003    -.0364711   -.0076785
       ivdem |   .1814579   .0503151     3.61   0.000     .0826384    .2802775
        time |  -.0007126   .0005768    -1.24   0.217    -.0018454    .0004201
       _cons |   -1.22033   .0365638   -33.38   0.000    -1.292142   -1.148518
------------------------------------------------------------------------------

.                  est store c3

.                  reghdfe v2lgotovst l1v2lgotovst l2v2lgotovst $d $x1,a(cowcode)cl
> uster(lid)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,364
Absorbing 1 HDFE group                            F(   6,    584) =     167.95
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9544
                                                  Adj R-squared   =     0.9522
                                                  Within R-sq.    =     0.6560
Number of clusters (lid)     =        585         Root MSE        =     0.2177

                                  (Std. err. adjusted for 585 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
  v2lgotovst | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
l1v2lgotovst |   .7805037   .0550284    14.18   0.000     .6724261    .8885814
l2v2lgotovst |  -.0094916   .0482801    -0.20   0.844    -.1043153    .0853322
   persparty |  -.1328246   .0389554    -3.41   0.001    -.2093343    -.056315
          ld |  -.0540801   .0185282    -2.92   0.004    -.0904701     -.01769
       ivdem |   .4640636   .1227634     3.78   0.000      .222952    .7051752
        time |   .0012513   .0009365     1.34   0.182     -.000588    .0030906
       _cons |  -1.129751   .0770556   -14.66   0.000    -1.281091   -.9784109
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       102           0         102     |
-----------------------------------------------------+

.                  est store d3

.                  xi:reghdfe v2lgotovst l1v2lgotovst l2v2lgotovst i.seat50*$d $x1,
> a(cowcode)cluster(lid)
i.seat50          _Iseat50_0-1        (naturally coded; _Iseat50_0 omitted)
i.seat50*pers~y   _IseaXpersp_#       (coded as above)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,190
Absorbing 1 HDFE group                            F(   8,    534) =     105.33
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9540
                                                  Adj R-squared   =     0.9517
                                                  Within R-sq.    =     0.6469
Number of clusters (lid)     =        535         Root MSE        =     0.2131

                                   (Std. err. adjusted for 535 clusters in lid)
-------------------------------------------------------------------------------
              |               Robust
   v2lgotovst | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
--------------+----------------------------------------------------------------
 l1v2lgotovst |   .7548006   .0593338    12.72   0.000     .6382443    .8713568
 l2v2lgotovst |   .0053969   .0505907     0.11   0.915    -.0939842    .1047781
   _Iseat50_1 |   .0324251   .0336753     0.96   0.336    -.0337272    .0985773
    persparty |  -.0886479   .0436938    -2.03   0.043    -.1744806   -.0028151
_IseaXpersp_1 |  -.1643199   .0791347    -2.08   0.038    -.3197734   -.0088663
           ld |  -.0507464   .0187843    -2.70   0.007    -.0876465   -.0138463
        ivdem |   .3912588   .1325011     2.95   0.003     .1309715    .6515461
         time |   .0011304   .0009457     1.20   0.233    -.0007274    .0029881
        _cons |  -1.068947   .0874645   -12.22   0.000    -1.240764   -.8971306
-------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       100           0         100     |
-----------------------------------------------------+

.                  est store s9

.                  interflex v2lgotovst $d v2paseat l1v2lgotovst l2v2lgotovst $x1,f
> e(cowcode)cluster(lid)nbin(4)
p value of Wald test: 0.0390

.                  mat list r(estBin)

r(estBin)[4,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1       17.15  -.06948686   .06420086  -.19531824   .05634452
r2        32.3  -.07513076   .05471053  -.18236142    .0320999
r3        45.1  -.13548684   .05423293  -.24178143  -.02919225
r4        56.6  -.25100967   .08101681   -.4097997  -.09221964

.                  reghdfe v2lgotovst l1v2lgotovst l2v2lgotovst $d $x1,a(cowcode)cl
> uster(lid)
(dropped 3 singleton observations)
(MWFE estimator converged in 1 iterations)

HDFE Linear regression                            Number of obs   =      2,364
Absorbing 1 HDFE group                            F(   6,    584) =     167.95
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9544
                                                  Adj R-squared   =     0.9522
                                                  Within R-sq.    =     0.6560
Number of clusters (lid)     =        585         Root MSE        =     0.2177

                                  (Std. err. adjusted for 585 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
  v2lgotovst | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
l1v2lgotovst |   .7805037   .0550284    14.18   0.000     .6724261    .8885814
l2v2lgotovst |  -.0094916   .0482801    -0.20   0.844    -.1043153    .0853322
   persparty |  -.1328246   .0389554    -3.41   0.001    -.2093343    -.056315
          ld |  -.0540801   .0185282    -2.92   0.004    -.0904701     -.01769
       ivdem |   .4640636   .1227634     3.78   0.000      .222952    .7051752
        time |   .0012513   .0009365     1.34   0.182     -.000588    .0030906
       _cons |  -1.129751   .0770556   -14.66   0.000    -1.281091   -.9784109
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |       102           0         102     |
-----------------------------------------------------+

.                  est store f3

.         
.                  erase .pdf

. 
.         **********************
.         *** Estimate plots ***
.         **********************
.                         gen e=.
(2,392 missing values generated)

.                         gen hi=.
(2,392 missing values generated)

.                         gen lo=.
(2,392 missing values generated)

.                         gen n =_n

.                         forval i =1/3   {               
  2.                                 est restore d`i'
  3.                                 nlcom _b[$d],post
  4.                                 mat e = e(b)
  5.                                 mat v = e(V)
  6.                                 mat list e
  7.                                 mat list v
  8.                                 local b=e[1,1] 
  9.                                 di `b'
 10.                                 local se = sqrt(v[1,1])
 11.                                 di `se'
 12.                                 replace e =`b' if n==`i'
 13.                                 replace hi = `b' + 1.96* `se' if n==`i'
 14.                                 replace lo = `b' - 1.96* `se'  if n==`i'
 15.                                 local i = `i'+1
 16.                         }
(results d1 are active now)

       _nl_1: _b[persparty]

------------------------------------------------------------------------------
v2xlg_legcon | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.1631885   .0425311    -3.84   0.000     -.246548    -.079829
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.16318847

symmetric v[1,1]
          _nl_1
_nl_1  .0018089
-.16318847
.04253114
(1 real change made)
(1 real change made)
(1 real change made)
(results d2 are active now)

       _nl_1: _b[persparty]

------------------------------------------------------------------------------
   v2x_jucon | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.0775784   .0310338    -2.50   0.012    -.1384035   -.0167533
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.07757843

symmetric v[1,1]
          _nl_1
_nl_1  .0009631
-.07757843
.03103377
(1 real change made)
(1 real change made)
(1 real change made)
(results d3 are active now)

       _nl_1: _b[persparty]

------------------------------------------------------------------------------
  v2lgotovst | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.1328246   .0389554    -3.41   0.001    -.2091757   -.0564736
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.13282465

symmetric v[1,1]
           _nl_1
_nl_1  .00151752
-.13282465
.03895535
(1 real change made)
(1 real change made)
(1 real change made)

.                  
.                         twoway (rspike hi lo n if n<=3,lcol(gs1)) (scatter e n if
>  n<=3,mcol(gs8)msym(O) ///
>                                 ytit({&beta}{sub:Party Personalist},size(large))y
> line(0,lcol(red)) ///
>                                 xlab(1 `""Legislative" "constraint""'   2  `""Jud
> icial" "constraint""' ///
>                                 3`""Bureaucratic" "oversight""' ) ///
>                                 xscale(range(0.8 3.2))xtit("")note(95% CI,size(vs
> mall)pos(8)ring(0))yscale(range(-.41 0.034))  legend(off) ///
>                                 xtit("Outcomes: institutional constraints on the 
> executive",height(6))saving(h1.gph,replace) ///
>                                 title(Party personalism reduces institutional con
> straints)ylab(-.4(.1)0))
(file h1.gph not found)
file h1.gph saved

.                         drop e n hi lo

.                 
.                         gen e=.
(2,392 missing values generated)

.                         gen hi=.
(2,392 missing values generated)

.                         gen lo=.
(2,392 missing values generated)

.                         gen n =_n

. 
.                         local var = "1 5 9"

.                         foreach i of local var   {      
  2.                                 est restore s`i'
  3.                                 nlcom _b[$d],post
  4.                                 mat e = e(b)
  5.                                 mat v = e(V)
  6.                                 mat list e
  7.                                 mat list v
  8.                                 local b=e[1,1] 
  9.                                 di `b'
 10.                                 local se = sqrt(v[1,1])
 11.                                 di `se'
 12.                                 replace e =`b' if n==`i'
 13.                                 replace hi = `b' + 1.96* `se' if n==`i'
 14.                                 replace lo = `b' - 1.96* `se'  if n==`i'
 15. 
.                                 est restore s`i'
 16.                                 nlcom _b[$d] + _b[_IseaXpersp_1],post
 17.                                 mat e = e(b)
 18.                                 mat v = e(V)
 19.                                 mat list e
 20.                                 mat list v
 21.                                 local b=e[1,1] 
 22.                                 di `b'
 23.                                 local se = sqrt(v[1,1])
 24.                                 di `se'
 25.                                 local j=`i'+1
 26.                                 replace e =`b' if n==`j'
 27.                                 replace hi = `b' + 1.96* `se' if n==`j'
 28.                                 replace lo = `b' - 1.96* `se'  if n==`j'
 29.                         }
(results s1 are active now)

       _nl_1: _b[persparty]

------------------------------------------------------------------------------
v2xlg_legcon | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.1360514   .0417791    -3.26   0.001     -.217937   -.0541659
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.13605143

symmetric v[1,1]
           _nl_1
_nl_1  .00174549
-.13605143
.04177911
(1 real change made)
(1 real change made)
(1 real change made)
(results s1 are active now)

       _nl_1: _b[persparty] + _b[_IseaXpersp_1]

------------------------------------------------------------------------------
v2xlg_legcon | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.2227356   .0834889    -2.67   0.008    -.3863709   -.0591003
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.22273562

symmetric v[1,1]
          _nl_1
_nl_1  .0069704
-.22273562
.08348895
(1 real change made)
(1 real change made)
(1 real change made)
(results s5 are active now)

       _nl_1: _b[persparty]

------------------------------------------------------------------------------
   v2x_jucon | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.0279127   .0312916    -0.89   0.372    -.0892431    .0334177
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.02791273

symmetric v[1,1]
           _nl_1
_nl_1  .00097916
-.02791273
.0312916
(1 real change made)
(1 real change made)
(1 real change made)
(results s5 are active now)

       _nl_1: _b[persparty] + _b[_IseaXpersp_1]

------------------------------------------------------------------------------
   v2x_jucon | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.2029035   .0601093    -3.38   0.001    -.3207155   -.0850915
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.20290349

symmetric v[1,1]
           _nl_1
_nl_1  .00361312
-.20290349
.06010926
(1 real change made)
(1 real change made)
(1 real change made)
(results s9 are active now)

       _nl_1: _b[persparty]

------------------------------------------------------------------------------
  v2lgotovst | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.0886479   .0436938    -2.03   0.042    -.1742861   -.0030096
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.08864785

symmetric v[1,1]
           _nl_1
_nl_1  .00190915
-.08864785
.04369377
(1 real change made)
(1 real change made)
(1 real change made)
(results s9 are active now)

       _nl_1: _b[persparty] + _b[_IseaXpersp_1]

------------------------------------------------------------------------------
  v2lgotovst | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -.2529677   .0723066    -3.50   0.000    -.3946861   -.1112493
------------------------------------------------------------------------------

symmetric e[1,1]
         _nl_1
y1  -.25296771

symmetric v[1,1]
           _nl_1
_nl_1  .00522825
-.25296771
.07230664
(1 real change made)
(1 real change made)
(1 real change made)

.                          
.                         twoway (rspike hi lo n if n==1 | n==5 | n==9  ,lcol(gs1))
>   ///
>                                 (rspike hi lo n if n==2 | n==6 | n==10  ,lcol(gs8
> ))  ///
>                                 (scatter e n if n==1 | n==5 | n==9  ,mcol(gs1)msy
> m(T)) ///
>                                 (scatter e n if n==2 | n==6 | n==10  ,mcol(gs8)ms
> ym(O) ///
>                                 ytit({&beta}{sub:Party Personalist},size(large))y
> line(0,lcol(red)) ///
>                                 xlab(1.5 `""Legislative" "constraint""' 5.5   `""
> Judicial" "constraint""' ///
>                                 9.5 `""Bureaucratic" "oversight""' ) ///
>                                 xscale(range(0.8 10.2))xtit("")note(95% CI,size(v
> small)pos(8)ring(0)) ///
>                                 legend(lab(3 "Ruling party seat share<50%")lab(4 
> "Ruling party seat share>=50%")pos(6)order(3 4)ring(0)col(2)) ///
>                                 xtit("Outcomes: institutional constraints on the 
> executive",height(6))yscale(range(-.41 0.034)) ///
>                                 title(By legislative majority) ///
>                                 ylab(-.4(.1)0) saving(h2.gph,replace))
(file h2.gph not found)
file h2.gph saved

.                         drop e n hi lo

.                         gr combine h1.gph h2.gph,xsize(9)tit(Party personalism an
> d institutional constraints on the executive)iscale(.7)

.                         gr export "$dir\golden\Ch5-Pers-Party-Exec-Constraint.pdf
> ",as(pdf)replace 
file
    C:\Users\jgw12\Dropbox\Research\PersPartyBook\Data\FKTW-reproduction\golden\C
    > h5-Pers-Party-Exec-Constraint.pdf saved as PDF format

.                         erase h1.gph

.                         erase h2.gph

. 
. 
.                         
. *************************************
. ****** Term limit analysis **********
. *************************************
. 
.                         *** Term limit changes ***
.                         use pers-use,clear

.                         tsset cowcode year

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit

.                         gen presvote = l.v2elvotlrg
(2,098 missing values generated)

.                         replace presvote = l.presvote if presvote==.
(990 real changes made)

.                         sort country year

.                         merge country year using termlimit
(you are using old merge syntax; see [D] merge for new syntax)

.                         tab _merge

     _merge |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |      2,364       97.20       97.20
          2 |         40        1.64       98.85
          3 |         28        1.15      100.00
------------+-----------------------------------
      Total |      2,432      100.00

.                         drop if _merge ==2
(40 observations deleted)

.                         replace termlimitattempt =1 if country=="Peru" & year==19
> 93
(0 real changes made)

.                         replace termlimitfail =0 if country=="Peru" & year==1993
(0 real changes made)

.                         replace termlimitattempt =1 if country=="Sri Lanka" & yea
> r==2010
(1 real change made)

.                         replace termlimitfail =0 if country=="Sri Lanka" & year==
> 2010
(1 real change made)

.                         recode termlimitfail termlimitattempt (.=0) if year<=2018
>  & year>=1993
(2,074 changes made to termlimitfail)
(2,074 changes made to termlimitattempt)

.                         gen election =v2xel_elecparl==1 |  v2xel_elecpres==1

.                         local var = "v2x_jucon v2xlg_legcon"

.                         foreach v of local var {
  2.                                 gen o`v'=l1`v' if year==min
  3.                                 egen i`v'=max(o`v'),by(lid)
  4.                         }
(1,815 missing values generated)
(55 missing values generated)
(1,817 missing values generated)
(71 missing values generated)

.                         xtset cowcode year

Panel variable: cowcode (unbalanced)
 Time variable: year, 1991 to 2020, but with gaps
         Delta: 1 unit

.                         gen l1election = l1.election
(132 missing values generated)

.                         gen l2election = l2.election
(251 missing values generated)

.                         gen f1election = f1.election
(132 missing values generated)

.                         gen f2election = f2.election
(251 missing values generated)

.                         gen ecycle = election==1 | f1election==1

.                         keep if year>=1993 & year<=2018
(289 observations deleted)

.                         tab country type

                      |         type
              country | parliam..  preside.. |     Total
----------------------+----------------------+----------
          Afghanistan |         0          4 |         4 
              Albania |        26          0 |        26 
            Argentina |         0         26 |        26 
              Armenia |         0          2 |         2 
            Australia |        26          0 |        26 
              Austria |        26          0 |        26 
           Azerbaijan |         0          1 |         1 
           Bangladesh |        20          0 |        20 
              Belgium |        26          0 |        26 
                Benin |         0         26 |        26 
              Bolivia |         0         26 |        26 
               Brazil |         0         26 |        26 
             Bulgaria |        26          0 |        26 
         Burkina Faso |         0          3 |         3 
              Burundi |         0          7 |         7 
               Canada |        26          0 |        26 
Central African Rep.. |         0         10 |        10 
                Chile |         0         26 |        26 
             Colombia |         0         26 |        26 
           Costa Rica |         0         26 |        26 
              Croatia |        26          0 |        26 
       Czech Republic |        23          0 |        23 
              Denmark |        26          0 |        26 
   Dominican Republic |         0         26 |        26 
              Ecuador |         0         26 |        26 
                Egypt |         0          1 |         1 
          El Salvador |         0         24 |        24 
              Estonia |        26          0 |        26 
Federal Republic of.. |         5          0 |         5 
              Finland |        18          8 |        26 
               France |         0         26 |        26 
               Gambia |         0          1 |         1 
              Georgia |         1         14 |        15 
              Germany |        26          0 |        26 
                Ghana |         0         18 |        18 
               Greece |        25          0 |        25 
            Guatemala |         0         23 |        23 
               Guinea |         0          8 |         8 
        Guinea Bissau |         0         13 |        13 
                Haiti |         0         17 |        17 
             Honduras |         0         26 |        26 
              Hungary |        26          0 |        26 
              Iceland |        26          0 |        26 
                India |        26          0 |        26 
            Indonesia |         0         19 |        19 
                 Iraq |         4          0 |         4 
              Ireland |        26          0 |        26 
               Israel |        26          0 |        26 
                Italy |        21          0 |        21 
          Ivory Coast |         0          7 |         7 
                Japan |        26          0 |        26 
                Kenya |         0         16 |        16 
               Kosovo |        10          0 |        10 
           Kyrgyzstan |         0          8 |         8 
               Latvia |        26          0 |        26 
              Lebanon |        13          0 |        13 
              Lesotho |        25          0 |        25 
              Liberia |         0         13 |        13 
            Lithuania |        26          0 |        26 
            Macedonia |        26          0 |        26 
           Madagascar |         0         21 |        21 
               Malawi |         0         24 |        24 
                 Mali |         0         25 |        25 
           Mauritania |         0          1 |         1 
            Mauritius |        26          0 |        26 
               Mexico |         0         18 |        18 
              Moldova |        26          0 |        26 
             Mongolia |         0         25 |        25 
                Nepal |        21          0 |        21 
          Netherlands |        26          0 |        26 
          New Zealand |        26          0 |        26 
            Nicaragua |         0         24 |        24 
                Niger |         0         20 |        20 
              Nigeria |         0         19 |        19 
               Norway |        26          0 |        26 
             Pakistan |        17          0 |        17 
               Panama |         0         26 |        26 
             Paraguay |         0         25 |        25 
                 Peru |         0         17 |        17 
          Philippines |         0         26 |        26 
               Poland |         0         26 |        26 
             Portugal |         0         26 |        26 
    Republic of Congo |         0          5 |         5 
              Romania |        25          0 |        25 
               Russia |         1          0 |         1 
              Senegal |         0         18 |        18 
               Serbia |        10          0 |        10 
    Serbia-Montenegro |         3          0 |         3 
         Sierra Leone |         0         21 |        21 
             Slovakia |        25          0 |        25 
             Slovenia |        26          0 |        26 
         South Africa |        24          0 |        24 
          South Korea |         0         26 |        26 
                Spain |        26          0 |        26 
            Sri Lanka |         0         19 |        19 
               Sweden |        26          0 |        26 
               Taiwan |         0         18 |        18 
             Thailand |        20          0 |        20 
              Tunisia |         4          3 |         7 
               Turkey |         0         24 |        24 
              Ukraine |         0         24 |        24 
       United Kingdom |        26          0 |        26 
        United States |         0         26 |        26 
              Uruguay |         0         26 |        26 
            Venezuela |         0         13 |        13 
               Zambia |         0         11 |        11 
----------------------+----------------------+----------
                Total |     1,043      1,060 |     2,103 

.                         keep if pres==1 | ccode==345
(1,025 observations deleted)

.                         
.                         * Bivariate descriptive *
.                         centile persparty,centile(50)

                                                          Binom. interp.   
    Variable |       Obs  Percentile    Centile        [95% conf. interval]
-------------+-------------------------------------------------------------
   persparty |     1,078         50    .6090533        .5948855    .6337572

.                         gen hip =persparty>=r(c_1) if persparty~=.              

.                         ttest termlimitatt,by(hip)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       0 |     539    .0185529    .0058177    .1350648    .0071248     .029981
       1 |     539    .0352505    .0079506    .1845835    .0196325    .0508684
---------+--------------------------------------------------------------------
Combined |   1,078    .0269017    .0049301    .1618712    .0172279    .0365755
---------+--------------------------------------------------------------------
    diff |           -.0166976    .0098517               -.0360284    .0026332
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.6949
H0: diff = 0                                     Degrees of freedom =     1076

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0452         Pr(|T| > |t|) = 0.0904          Pr(T > t) = 0.9548

.                                 
.                          * Variables *
.                         alpha presvote v2paseat v2pavote,std item gen(prespopular
> )

Test scale = mean(standardized items)

                                                            Average
                             Item-test     Item-rest       interitem
Item         |  Obs  Sign   correlation   correlation     correlation     alpha
-------------+-----------------------------------------------------------------
presvote     |  836    +       0.7670        0.4359          0.9176      0.9570
v2paseatsh~e |  999    +       0.9114        0.7284          0.4683      0.6379
v2pavote     |  838    +       0.9420        0.8450          0.4022      0.5737
-------------+-----------------------------------------------------------------
Test scale   |                                               0.6108      0.8248
-------------------------------------------------------------------------------

.                         local var = "i_popul persparty ivdem iv2x_jucon iv2xlg_le
> gcon ld"

.                         foreach v of local var {
  2.                                 qui sum `v'
  3.                                 replace `v' = (`v' - abs(r(mean)))/r(sd)
  4.                          
.                         }
(1,029 real changes made)
(1,078 real changes made)
(1,078 real changes made)
(1,062 real changes made)
(1,048 real changes made)
(1,078 real changes made)

.                         replace ivdem  =ivdem*10 /* rescale to plot estimate */
(1,078 real changes made)

.                         sum i_pop persparty

    Variable |        Obs        Mean    Std. dev.       Min        Max
-------------+---------------------------------------------------------
  i_populism |      1,029    3.78e-09           1  -1.676303   2.196776
   persparty |      1,078    3.02e-09           1  -2.770436   1.956474

.                         replace l1v2x_jucon=l1v2x_jucon*2
(1,078 real changes made)

.                         replace l1v2xlg_legcon=l1v2xlg_legcon*2
(1,073 real changes made)

.                         gen Election = election==1 | f1election==1

.                         
.                         twoway (lpoly termlimitattempt persparty,legend(off)bw(.5
> ))  (lowess termlimitattempt persparty) 

.                         probit termlimitattempt ld persparty,cluster(lid)

Iteration 0:  Log pseudolikelihood = -133.45782  
Iteration 1:  Log pseudolikelihood = -128.60011  
Iteration 2:  Log pseudolikelihood = -128.43099  
Iteration 3:  Log pseudolikelihood = -128.43051  
Iteration 4:  Log pseudolikelihood = -128.43051  

Probit regression                                       Number of obs =  1,078
                                                        Wald chi2(2)  =   9.52
                                                        Prob > chi2   = 0.0086
Log pseudolikelihood = -128.43051                       Pseudo R2     = 0.0377

                                      (Std. err. adjusted for 238 clusters in lid)
----------------------------------------------------------------------------------
                 |               Robust
termlimitattempt | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
-----------------+----------------------------------------------------------------
              ld |   .2599467   .0842889     3.08   0.002     .0947435    .4251499
       persparty |   .1891626   .0973528     1.94   0.052    -.0016453    .3799706
           _cons |   -1.99681   .0920627   -21.69   0.000    -2.177249    -1.81637
----------------------------------------------------------------------------------

. 
.                          * Reported estimates *
.                         reghdfe termlimitattempt ld persparty,a(cowcode year)clus
> ter(lid)
(dropped 4 singleton observations)
(MWFE estimator converged in 6 iterations)

HDFE Linear regression                            Number of obs   =      1,074
Absorbing 2 HDFE groups                           F(   2,    233) =       3.74
Statistics robust to heteroskedasticity           Prob > F        =     0.0251
                                                  R-squared       =     0.0870
                                                  Adj R-squared   =     0.0104
                                                  Within R-sq.    =     0.0044
Number of clusters (lid)     =        234         Root MSE        =     0.1613

                                  (Std. err. adjusted for 234 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
termlimita~t | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
          ld |   .0116118   .0076418     1.52   0.130     -.003444    .0266677
   persparty |   .0150024   .0066371     2.26   0.025     .0019261    .0280787
       _cons |   .0269348   .0045064     5.98   0.000     .0180562    .0358134
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        57           0          57     |
        year |        26           1          25     |
-----------------------------------------------------+

.                         est store term1

.                         gen s = e(sample)==1

.                         reghdfe termlimitattempt ivdem Election ld persparty,a(co
> wcode year)cluster(lid)
(dropped 4 singleton observations)
(MWFE estimator converged in 6 iterations)

HDFE Linear regression                            Number of obs   =      1,074
Absorbing 2 HDFE groups                           F(   4,    233) =       3.62
Statistics robust to heteroskedasticity           Prob > F        =     0.0070
                                                  R-squared       =     0.0940
                                                  Adj R-squared   =     0.0161
                                                  Within R-sq.    =     0.0122
Number of clusters (lid)     =        234         Root MSE        =     0.1609

                                  (Std. err. adjusted for 234 clusters in lid)
------------------------------------------------------------------------------
             |               Robust
termlimita~t | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       ivdem |  -.0019167   .0017145    -1.12   0.265    -.0052947    .0014612
    Election |   .0275812   .0091745     3.01   0.003     .0095056    .0456568
          ld |   .0089127   .0077798     1.15   0.253    -.0064151    .0242405
   persparty |   .0165945   .0069716     2.38   0.018      .002859    .0303299
       _cons |   .0115673   .0050739     2.28   0.024     .0015707    .0215639
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        57           0          57     |
        year |        26           1          25     |
-----------------------------------------------------+

.                         est store term2

.                         reghdfe termlimitattempt ivdem prespopular Election iv2x_
> jucon iv2xlg_legcon ld persparty,a(cowcode year)cluster(lid)
(dropped 3 singleton observations)
(MWFE estimator converged in 6 iterations)

HDFE Linear regression                            Number of obs   =      1,022
Absorbing 2 HDFE groups                           F(   7,    221) =       2.51
Statistics robust to heteroskedasticity           Prob > F        =     0.0166
                                                  R-squared       =     0.1007
                                                  Adj R-squared   =     0.0179
                                                  Within R-sq.    =     0.0173
Number of clusters (lid)     =        222         Root MSE        =     0.1646

                                   (Std. err. adjusted for 222 clusters in lid)
-------------------------------------------------------------------------------
              |               Robust
termlimitat~t | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
--------------+----------------------------------------------------------------
        ivdem |  -.0002089   .0023045    -0.09   0.928    -.0047506    .0043327
  prespopular |   .0086854   .0094574     0.92   0.359    -.0099528    .0273235
     Election |   .0308225   .0098051     3.14   0.002      .011499     .050146
   iv2x_jucon |  -.0233007   .0172358    -1.35   0.178    -.0572682    .0106667
iv2xlg_legcon |  -.0094028   .0124054    -0.76   0.449    -.0338508    .0150451
           ld |   .0105083   .0108364     0.97   0.333    -.0108476    .0318642
    persparty |   .0197387   .0074876     2.64   0.009     .0049825     .034495
        _cons |   .0117021   .0055703     2.10   0.037     .0007245    .0226797
-------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        55           0          55     |
        year |        26           1          25     |
-----------------------------------------------------+

.                         est store term3

.                         reghdfe termlimitattempt i_populism ivdem prespopular Ele
> ction iv2x_jucon iv2xlg_legcon ld persparty,a(cowcode year)cluster(lid)
(dropped 3 singleton observations)
(MWFE estimator converged in 6 iterations)

HDFE Linear regression                            Number of obs   =        990
Absorbing 2 HDFE groups                           F(   8,    214) =       2.59
Statistics robust to heteroskedasticity           Prob > F        =     0.0101
                                                  R-squared       =     0.1044
                                                  Adj R-squared   =     0.0180
                                                  Within R-sq.    =     0.0205
Number of clusters (lid)     =        215         Root MSE        =     0.1672

                                   (Std. err. adjusted for 215 clusters in lid)
-------------------------------------------------------------------------------
              |               Robust
termlimitat~t | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
--------------+----------------------------------------------------------------
   i_populism |   .0102016   .0069928     1.46   0.146     -.003582    .0239852
        ivdem |  -.0019182   .0026488    -0.72   0.470    -.0071393     .003303
  prespopular |   .0070022   .0094276     0.74   0.458    -.0115805     .025585
     Election |   .0315211   .0100656     3.13   0.002     .0116806    .0513615
   iv2x_jucon |  -.0264779   .0168266    -1.57   0.117    -.0596449    .0066892
iv2xlg_legcon |  -.0072386   .0129779    -0.56   0.578    -.0328194    .0183422
           ld |   .0125467   .0110724     1.13   0.258    -.0092782    .0343716
    persparty |     .01833   .0075847     2.42   0.017     .0033797    .0332803
        _cons |   .0146201   .0059426     2.46   0.015     .0029065    .0263337
-------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
     cowcode |        55           0          55     |
        year |        26           1          25     |
-----------------------------------------------------+

.                         est store term4

.                         
.                          * Plot estimates *
.                         label var iv2x_jucon `""Judicial  " "constraint""'

.                         label var iv2xlg_legcon `""Legislative" "constraint ""'

.                         label var persparty `""{bf:Party}     " "{bf:personalism}
> ""'

.                         label var ivdem `""Initial    " "democracy" "level     ""
> '

.                         label var ld `""Democracy" "age     ""'

.                         label var Election "Election year"

.                         label var i_populism `""Party   " "Populism""'

.                         label var prespopular `""Leader  " "Popularity""'

. 
.                         coefplot (term1, msymbol(d))(term2, msymbol(P)) (term3, m
> symbol(T)) (term4, msymbol(Oh)), order(persparty)  ///
>                                 drop(_cons) xline(0) msymbol(d) mfcolor(white) gr
> id(glcolor(gs15)) ///
>                                 levels(95 90) legend(lab(3 "Baseline")lab(6 "+ el
> ection + initial democracy level") ///
>                                 lab(9 "+ legislative/judicial constraint")  lab(1
> 2 "+ populism")order(3 6 9 12) ///
>                                 size(small) pos(6) col(2) ring(1)) xsize(2) ysize
> (2) xlab(-.06(.02).04)  ///
>                                 xtitle("        Coefficient estimate", size(small
> ))  ///
>                                 ciopts(lwidth(thin)) aspectratio(1.1) scale(.75) 
> title(Attempts to adjust term limits, size(medium) height(2)) ///
>                                 note("Presidential systems only",pos(6)size(vsmal
> l))
(note:  named style P not found in class symbol, default attributes used)

.                         gr export "$dir\golden\Ch5-Pers-Party-Term-Limit-Attempts
> .pdf",as(pdf)replace 
file
    C:\Users\jgw12\Dropbox\Research\PersPartyBook\Data\FKTW-reproduction\golden\C
    > h5-Pers-Party-Term-Limit-Attempts.pdf saved as PDF format

.         
.                         * Check with non-linear link *
.                         xthybrid termlimitattempt persparty ld  if s==1,cluster(c
> owcode)fam(binomial)link(logit)p


Hybrid model. Family: binomial. Link: logit.

+-----------------------------------+
|             Variable |   model    |
|----------------------+------------|
| termlimitattempt     |            |
|         W__persparty |     0.5935 |
|                      |     0.0625 |
|                W__ld |     1.0985 |
|                      |     0.0251 |
|         B__persparty |     0.2547 |
|                      |     0.4542 |
|                B__ld |     0.4816 |
|                      |     0.1097 |
|                _cons |    -4.0701 |
|                      |     0.0000 |
|----------------------+------------|
|   var(_cons[cowcode])|            |
|                _cons |     0.5001 |
|                      |     0.2904 |
|----------------------+------------|
| Statistics           |            |
|                   ll |  -126.4881 |
|                 chi2 |     9.6310 |
|                    p |     0.0471 |
|                  aic |   264.9762 |
|                  bic |   294.8510 |
+-----------------------------------+
                          Legend: b/p
Level 1: 1074 units. Level 2: 57 units.

.                         xthybrid termlimitattempt persparty xperiod* Election ld 
> ivdem  if s==1,cluster(cowcode)fam(binomial)link(logit)p


Hybrid model. Family: binomial. Link: logit.

+-----------------------------------+
|             Variable |   model    |
|----------------------+------------|
| termlimitattempt     |            |
|         W__persparty |     0.6965 |
|                      |     0.0435 |
|          W__xperiod1 |     0.1728 |
|                      |     0.9134 |
|          W__xperiod2 |     0.9365 |
|                      |     0.4165 |
|          W__xperiod3 |    -0.3260 |
|                      |     0.7725 |
|          W__xperiod4 |     1.1397 |
|                      |     0.1399 |
|          W__xperiod5 |     0.9075 |
|                      |     0.2142 |
|          W__xperiod6 |  (omitted) |
|                      |            |
|          W__Election |     1.2938 |
|                      |     0.0088 |
|                W__ld |     1.1164 |
|                      |     0.2031 |
|             W__ivdem |    -0.0623 |
|                      |     0.3933 |
|         B__persparty |     0.0820 |
|                      |     0.8263 |
|          B__xperiod1 |     6.2155 |
|                      |     0.4558 |
|          B__xperiod2 |     2.6371 |
|                      |     0.7552 |
|          B__xperiod3 |     9.0035 |
|                      |     0.3125 |
|          B__xperiod4 |     4.6546 |
|                      |     0.6097 |
|          B__xperiod5 |     4.4213 |
|                      |     0.7008 |
|          B__xperiod6 |  (omitted) |
|                      |            |
|          B__Election |     0.1741 |
|                      |     0.9103 |
|                B__ld |     0.5931 |
|                      |     0.1824 |
|             B__ivdem |    -0.0338 |
|                      |     0.3274 |
|                _cons |    -9.2309 |
|                      |     0.2416 |
|----------------------+------------|
|   var(_cons[cowcode])|            |
|                _cons |     0.5207 |
|                      |     0.3087 |
|----------------------+------------|
| Statistics           |            |
|                   ll |  -116.9398 |
|                 chi2 |    22.5119 |
|                    p |     0.2100 |
|                  aic |   273.8796 |
|                  bic |   373.4625 |
+-----------------------------------+
                          Legend: b/p
Level 1: 1074 units. Level 2: 57 units.

.                         xthybrid termlimitattempt persparty xperiod* ivdem prespo
> pular Election iv2x_jucon iv2xlg_legcon ld ///
>                                 if s==1,cluster(cowcode)fam(binomial)link(logit)p


Hybrid model. Family: binomial. Link: logit.

+-----------------------------------+
|             Variable |   model    |
|----------------------+------------|
| termlimitattempt     |            |
|         W__persparty |     0.8763 |
|                      |     0.0254 |
|          W__xperiod1 |     0.7988 |
|                      |     0.6419 |
|          W__xperiod2 |     1.1433 |
|                      |     0.3764 |
|          W__xperiod3 |    -0.4393 |
|                      |     0.7239 |
|          W__xperiod4 |     1.1318 |
|                      |     0.1725 |
|          W__xperiod5 |     0.8051 |
|                      |     0.2814 |
|          W__xperiod6 |  (omitted) |
|                      |            |
|             W__ivdem |     0.0235 |
|                      |     0.7980 |
|       W__prespopular |     0.2205 |
|                      |     0.4982 |
|          W__Election |     1.3043 |
|                      |     0.0080 |
|        W__iv2x_jucon |    -0.8788 |
|                      |     0.3186 |
|     W__iv2xlg_legcon |    -0.2698 |
|                      |     0.6736 |
|                W__ld |     1.2094 |
|                      |     0.2291 |
|         B__persparty |    -0.0865 |
|                      |     0.8281 |
|          B__xperiod1 |    -3.3320 |
|                      |     0.6924 |
|          B__xperiod2 |     0.2976 |
|                      |     0.9666 |
|          B__xperiod3 |     5.1403 |
|                      |     0.4477 |
|          B__xperiod4 |    -0.5895 |
|                      |     0.9219 |
|          B__xperiod5 |     0.0879 |
|                      |     0.9927 |
|          B__xperiod6 |  (omitted) |
|                      |            |
|             B__ivdem |     0.0165 |
|                      |     0.6979 |
|       B__prespopular |     0.0946 |
|                      |     0.8246 |
|          B__Election |    -1.3174 |
|                      |     0.3608 |
|        B__iv2x_jucon |    -0.6497 |
|                      |     0.0599 |
|     B__iv2xlg_legcon |    -0.3149 |
|                      |     0.3087 |
|                B__ld |     0.7881 |
|                      |     0.0796 |
|                _cons |    -4.3828 |
|                      |     0.4605 |
|----------------------+------------|
|   var(_cons[cowcode])|            |
|                _cons |     0.0284 |
|                      |     0.9465 |
|----------------------+------------|
| Statistics           |            |
|                   ll |  -109.9582 |
|                 chi2 |    32.9850 |
|                    p |     0.1044 |
|                  aic |   271.9165 |
|                  bic |   400.0839 |
+-----------------------------------+
                          Legend: b/p
Level 1: 1022 units. Level 2: 55 units.

.                         xthybrid termlimitattempt persparty xperiod* i_popul ivde
> m prespopular Election iv2x_jucon iv2xlg_legcon ld ///
>                                 if s==1,cluster(cowcode)fam(binomial)link(logit)p


Hybrid model. Family: binomial. Link: logit.

+-----------------------------------+
|             Variable |   model    |
|----------------------+------------|
| termlimitattempt     |            |
|         W__persparty |     0.7842 |
|                      |     0.0491 |
|          W__xperiod1 |     0.8832 |
|                      |     0.6102 |
|          W__xperiod2 |     1.2846 |
|                      |     0.3259 |
|          W__xperiod3 |    -0.3409 |
|                      |     0.7875 |
|          W__xperiod4 |     1.0962 |
|                      |     0.1872 |
|          W__xperiod5 |     0.6923 |
|                      |     0.3562 |
|          W__xperiod6 |  (omitted) |
|                      |            |
|        W__i_populism |     0.1789 |
|                      |     0.5644 |
|             W__ivdem |    -0.0194 |
|                      |     0.8447 |
|       W__prespopular |     0.1952 |
|                      |     0.5674 |
|          W__Election |     1.2995 |
|                      |     0.0084 |
|        W__iv2x_jucon |    -0.7850 |
|                      |     0.3890 |
|     W__iv2xlg_legcon |    -0.2513 |
|                      |     0.7074 |
|                W__ld |     1.3832 |
|                      |     0.1907 |
|         B__persparty |    -0.2138 |
|                      |     0.6452 |
|          B__xperiod1 |    -2.2189 |
|                      |     0.7907 |
|          B__xperiod2 |     6.8902 |
|                      |     0.4048 |
|          B__xperiod3 |     9.2173 |
|                      |     0.1722 |
|          B__xperiod4 |     4.0138 |
|                      |     0.5559 |
|          B__xperiod5 |     7.9548 |
|                      |     0.4518 |
|          B__xperiod6 |  (omitted) |
|                      |            |
|        B__i_populism |     0.3435 |
|                      |     0.4587 |
|             B__ivdem |    -0.0058 |
|                      |     0.9104 |
|       B__prespopular |    -0.1181 |
|                      |     0.7898 |
|          B__Election |    -1.2888 |
|                      |     0.3766 |
|        B__iv2x_jucon |    -0.5150 |
|                      |     0.1482 |
|     B__iv2xlg_legcon |    -0.2801 |
|                      |     0.3954 |
|                B__ld |     0.8618 |
|                      |     0.0699 |
|                _cons |    -9.0620 |
|                      |     0.1642 |
|----------------------+------------|
|   var(_cons[cowcode])|            |
|                _cons |     0.0236 |
|                      |     0.9544 |
|----------------------+------------|
| Statistics           |            |
|                   ll |  -108.0590 |
|                 chi2 |    33.1434 |
|                    p |     0.1580 |
|                  aic |   272.1181 |
|                  bic |   409.2538 |
+-----------------------------------+
                          Legend: b/p
Level 1: 990 units. Level 2: 55 units.

.         
.         
.         ** By legislative majority **
.         ttest  termlimitattempt  if v2paseat~=. & l1election==0,by(create)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       0 |     431    .0324826    .0085491     .177484    .0156794    .0492858
       1 |     240    .0416667    .0129257    .2002439    .0162038    .0671295
---------+--------------------------------------------------------------------
Combined |     671    .0357675    .0071746    .1858485    .0216801    .0498549
---------+--------------------------------------------------------------------
    diff |           -.0091841    .0149754               -.0385885    .0202204
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.6133
H0: diff = 0                                     Degrees of freedom =      669

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2700         Pr(|T| > |t|) = 0.5399          Pr(T > t) = 0.7300

.         ttest  termlimitattempt  if v2paseats< 50 & v2paseat~=. & l1election==0,b
> y(create)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       0 |     280    .0285714     .009974    .1668969    .0089376    .0482053
       1 |     181    .0276243    .0122159    .1643485    .0035195    .0517292
---------+--------------------------------------------------------------------
Combined |     461    .0281996    .0077185    .1657224    .0130317    .0433674
---------+--------------------------------------------------------------------
    diff |            .0009471    .0158228               -.0301471    .0320413
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.0599
H0: diff = 0                                     Degrees of freedom =      459

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.5239         Pr(|T| > |t|) = 0.9523          Pr(T > t) = 0.4761

.         ttest  termlimitattempt  if v2paseats>=50 & v2paseat~=. & l1election==0,b
> y(create)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       0 |     151    .0397351    .0159491    .1959862    .0082211    .0712491
       1 |      59    .0847458    .0365692    .2808936    .0115445     .157947
---------+--------------------------------------------------------------------
Combined |     210     .052381     .015411    .2233264        .022    .0827619
---------+--------------------------------------------------------------------
    diff |           -.0450107    .0342278               -.1124885    .0224672
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.3150
H0: diff = 0                                     Degrees of freedom =      208

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0950         Pr(|T| > |t|) = 0.1899          Pr(T > t) = 0.9050

.         
.         interflex termlimitattempt persparty v2paseat ld ivdem if l1election==0,f
> e(cowcode year)cluster(lid)nbin(2)cutoffs(49.9)
p value of Wald test: 0.0672

.         mat list r(estBin)

r(estBin)[2,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        32.7   .01139787    .0118215  -.01177184   .03456759
r2        56.3   .06334088   .02019603   .02375739   .10292437

.         interflex termlimitattempt create v2paseat ld ivdem if l1election==0,fe(c
> owcode year)cluster(lid)nbin(2)cutoffs(49.9)
p value of Wald test: 0.1131

.         mat list r(estBin)

r(estBin)[2,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        32.7   .01424547   .02179582  -.02847356    .0569645
r2        56.3   .10558319   .04374839   .01983792   .19132846

.         interflex termlimitattempt persparty v2paseat ld ivdem if l1election==0,f
> e(cowcode year)cluster(lid)nbin(2)cutoffs(49.9)
p value of Wald test: 0.0672

.         mat list r(estBin)

r(estBin)[2,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        32.7   .01139787    .0118215  -.01177184   .03456759
r2        56.3   .06334088   .02019603   .02375739   .10292437

.         interflex termlimitattempt create v2paseat ld ivdem if l1election==0,fe(c
> owcode year)cluster(lid)nbin(2)cutoffs(49.9)
p value of Wald test: 0.1131

.         mat list r(estBin)

r(estBin)[2,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        32.7   .01424547   .02179582  -.02847356    .0569645
r2        56.3   .10558319   .04374839   .01983792   .19132846

.         interflex termlimitattempt persparty v2paseat i_populism ivdem prespopula
> r Election iv2x_jucon iv2xlg_legcon ld ///
>                 if l1election==0,fe(cowcode year)cluster(lid)cut(49.9)
p value of Wald test: 0.0433

.         mat list r(estBin)

r(estBin)[2,5]
           x0   bin_marg     bin_se   bin_CI_l   bin_CI_u
r1       32.3  .00984606  .01179591  -.0132735  .03296563
r2       56.3  .06659705  .02135862  .02473493  .10845917

.         interflex termlimitattempt create v2paseat i_populism ivdem prespopular E
> lection iv2x_jucon iv2xlg_legcon ld ///
>                 if l1election==0,fe(cowcode year)cluster(lid)cut(49.9)
p value of Wald test: 0.1129

.         mat list r(estBin)

r(estBin)[2,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        32.3   .01277569   .02312857  -.03255547   .05810685
r2        56.3   .07754927   .03870961   .00167983   .15341871

. 
.         * Result reported in online Appendix *
.         gen e=.
(1,078 missing values generated)

.         gen hi=.
(1,078 missing values generated)

.         gen lo=.
(1,078 missing values generated)

.         gen n =_n

.         interflex termlimitattempt persparty v2paseat ld ivdem if l1election==0,f
> e(cowcode year)cluster(lid)cut(49.9)
p value of Wald test: 0.0672

.         mat list r(estBin)

r(estBin)[2,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        32.7   .01139787    .0118215  -.01177184   .03456759
r2        56.3   .06334088   .02019603   .02375739   .10292437

. 
.         replace e=r(estBin)[1,2] if n==1
(1 real change made)

.         replace e=r(estBin)[2,2] if n==2
(1 real change made)

.         replace hi=r(estBin)[1,5] if n==1
(1 real change made)

.         replace hi=r(estBin)[2,5] if n==2
(1 real change made)

.         replace lo=r(estBin)[1,4] if n==1
(1 real change made)

.         replace lo=r(estBin)[2,4] if n==2
(1 real change made)

.         twoway (rspike lo hi n if n<=2,legend(off)ytit(Marginal effect of ruling 
> party personalism) xtit(Ruling party legislative seat share))  ///
>                 (scatter e n if n<=2,xlab(1 "<50 %" 2 ">=50%")yline(0)xscale(rang
> e(0.8 2.2))saving(h1.gph,replace)tit(Ruling party personalism) ///
>                 yscale(range(-0.04 .2))ylab(0(.05).2)) 
(file h1.gph not found)
file h1.gph saved

.                 
.                 
.         interflex termlimitattempt create v2paseat ld ivdem if l1election==0,fe(c
> owcode year)cluster(lid)cut(49.9)
p value of Wald test: 0.1131

.         mat list r(estBin)

r(estBin)[2,5]
            x0    bin_marg      bin_se    bin_CI_l    bin_CI_u
r1        32.7   .01424547   .02179582  -.02847356    .0569645
r2        56.3   .10558319   .04374839   .01983792   .19132846

. 
.         replace e=r(estBin)[1,2] if n==1
(1 real change made)

.         replace e=r(estBin)[2,2] if n==2
(1 real change made)

.         replace hi=r(estBin)[1,5] if n==1
(1 real change made)

.         replace hi=r(estBin)[2,5] if n==2
(1 real change made)

.         replace lo=r(estBin)[1,4] if n==1
(1 real change made)

.         replace lo=r(estBin)[2,4] if n==2
(1 real change made)

.         twoway (rspike lo hi n if n<=2,legend(off)ytit(Marginal effect of create-
> party) xtit(Ruling party legislative seat share))  ///
>                 (scatter e n if n<=2,xlab(1 "<50 %" 2 ">=50%")yline(0)xscale(rang
> e(0.8 2.2))saving(h2.gph,replace)tit(Create party) ///
>                 yscale(range(-0.04 .2))ylab(0(.05).2)) 
(file h2.gph not found)
file h2.gph saved

.                 
.         gr combine h1.gph h2.gph,tit(Party personalism increases term limit attem
> pts when ruling party has legislative majority)xsize(8)

.         gr export "$dir\golden\T-Term-Limit-Attempts-LegSeatShare.pdf",as(pdf)rep
> lace 
file
    C:\Users\jgw12\Dropbox\Research\PersPartyBook\Data\FKTW-reproduction\golden\T
    > -Term-Limit-Attempts-LegSeatShare.pdf saved as PDF format

.         
.         erase .pdf

.         erase h1.gph

.         erase h2.gph

. 
.         ************  THE END ***********
.                         
.  log close
      name:  <unnamed>
       log:  C:\Users\jgw12\Dropbox\Research\PersPartyBook\Data\FKTW-reproduction\C
> h5.log
  log type:  text
 closed on:  26 Jul 2023, 17:04:15
-----------------------------------------------------------------------------------
